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We analyze all publicly available long-term optical observations of the gravita¬ 
tionally lensed quasar PG1115+080 for the purpose of estimating time delays between 
its four components. In particular, the light curves of PGl 115+080 components 
obtained in 2001-2006 at Maidanak observatory (Uzbekistan) [60] are considered. We 
hnd that the linear trend is observed in 2006 in light curves of all four components 
with fast variations only in the A 1 and G components that can be due to microlensing 
and observational errors. Application of the MGCP method [36] to the photometric 
data obtained in 2004-2005 gives values of time delays tbc = 22 ^( 3 , tac = 12^^ 
and tba = IOI 3 days, which are in agreement with the results received earlier by 
Schechter [48] and Barkana [2] for 1995-1996 light curves with two different methods 
of statistic analysis. However, our estimates of tba and tbc differ from the values 
received by the group of Vakulik based on the same Maidanak data [61]. The ratio 
tac/tba is equal to ~ 1.2 that is close to the value, received by Barkana (~ 1.13) 
and predicted by lens models (~ 1-4), unlike the values received by Schechter (~ 0.7) 
and Vakulik (~ 2.7). 
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1. INTRODUCTION 

According to General relativity, beams of light deviate from the straight line under the 
influence of the gravitational held of a massive object [11]. If a galaxy is located close to the 
line of sight between a far quasar and an observer, the light from the quasar passes different 
ways that leads to formation of several images of the quasar. The internal variability of a 
quasar manifests itself in light curves of all these images, but with time shifts. The time 
delay between manifestation of the internal variability of a quasar in a couple of its images 
is caused by geometrical difference of ways, light pass for each image, and difference of the 
gravitational potential in points corresponding to positions of quasar images. 

Even before the discovery of the hrst gravitational lens, the double quasar QSO 0957+561, 
Refsdal showed that the time delay between manifestations of hux changes of a source (a 
supernova or a quasar) which is located behind a remote lensing galaxy and close to the line 
of sight allows determining the Hubble constant Hq for certain model of the gravitational 
potential of a lensing galaxy [43, 44]. And vice versa, it is possible to investigate models of a 
lensing galaxy mass distribution using measured time delays and an estimated value of the 
Hubble constant, obtained with other methods. 

A number of discovered gravitational lenses is increasing and every new one can become a 
candidate for Hubble constant estimation that allows minimizing systematic errors. By now, 
delays between components are measured for more than 20 lenses. Difficulties of time delay 
measurement in gravitational lens systems are concerned with the fact that the majority 
of gravitationally lensed quasars (GLQ) has rather small internal variability and there are 
manifestations of the micro-lensing effect due to separate stars of a lensing galaxy passing 
near ways of light forming images of a quasar. Therefore long-term homogeneous observations 
of GLQs with high resolution and robust time delay measurement methods are needed. 

In this work, we analyze light curves of the components of the gravitationally lensed 
quasar PG 1115+080 to determine time delays. This is the second GLQ discovered and the 
hrst quadruple one. During hrst observations, three components A, B and G with relatively 
large angle separation (about 2") were discovered [64]. Later the brightest component A was 
resolved into two components Al and A2 spaced 0."48 apart [16], that conhrmed assumptions 
by Young et al [65]. The system conhguration corresponds to the source position near a 
macrocaustics fold. The redshift of the source is zg = 1.722 and the redshift of the lensing 
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galaxy is zq = 0.31 [7, 17, 59]. Thanks to its well-studied geometry [10, 27], evident source 
variability and theoretically estimated time delays not exceeding dozens of days [30, 63], this 
system became the second GLS that was used to estimate the Hubble constant [ 22 , 48]. 

Long-term observation sets are needed to determine time delays. Such time series for 
PG1115+080 were hrst published by Schechter [48] in 1997. The quasar was observed in 
V band from November 1995 to June 1996. The Press method [42], used to estimate time 
delays, showed that flux variations of B component follow flux variations of G component 
with the delay of 23.7 ± 3.4 days, and flux variations of A1 and A2 follow flux variations of 
C with the delay of 9.4 ± 3.4 days. It is notable that time delay between Al and A 2 is less 
than a day. X-ray observations allowed estimating taiA 2 = 0.149 ± 0.004 days [ 6 ]. Therefore 
the light curves of Ai and A 2 are usually combined for time delay estimation. 

The same data were analyzed by Barkana [2] with different statistical approach. Using 
analytical representation of the light curves and taking into account correlated errors of 
photometry, he obtained tcb = 25.0^3 § days, which is in agreement with Schechter’s estimate 
within the accuracy of the analysis, and tca = 13.412 2 days, which is much larger then 
Schechter’s value. Another method of time delay estimation was applied to the same data 
set by Pelt et al. The time delays calculated with the minimum dispersion method are: 
Tab = 15.5 ± 1.8, tca = 10.3 ± 1.9, tcb = 25.8 ± 2.4 days [40]. 

New estimates of time delays in PG1115+080 were published in 2009 based on the analysis 
of images obtained with the 1.5-m telescope of the Maidanak Observatory (Uzbekistan) 
[61] in 2004-2006. The authors applied their time delay estimation method and reported 
Tcb = 16 . 41^5 and tac = 12.Ol^^ days. Those values substantially differ from Schechter’s 
and Barkana’s results and may suggest higher value of the Hubble constant. 

We analysed the same Maidanak data sets with our modihcation of the Gaskell and Sparke 
method of the cross-correlation analysis [14, 15, 36] and obtained results supporting ’old’ 
values of Schechter and Barkana [1]. 

Another attempt to measure time delays in PG1115+080 based on observations in 1995- 
1996 was undertaken in 2011 by Eulaers and Magain using two fundamentally different 
techniques [12]. Modihed numerical model £t technique [5] resulted in tqb = 23.8^3 0 and 
Tca = 11-7 ± 2.2, and minimum dispersion technique [39] resulted in tcb = 17.9 ± 6.9 and 
Tca = 7.6 ± 3.9 days. Because of such difference they concluded that estimation of time 
delays based on those data is impossible. 
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Recently the Maidanak observations were re-analysed by Rathna Kumar et al. using the 
difference-smoothing technique that resulted in tcb = 13.2 ± 9.0 and tca = 18.3 ± 4.4. 

As the discrepancy of time delays obtained with different techniques has significant 
impact on the Hubble constant estimation, we aimed to more thoroughly inspect available 
observational data and apply our technique of time delay estimation. In the following 
sections, we analyze available time series for PG1115+080, describe our modification of the 
Gaskell and Sparke method of the cross-correlation analysis [14, 15, 36] and discuss results 
of its application to publicly available observational data sets. 

2. OBSERVATIONS AND REDUGTION 

There are three available long-term monitoring data sets for PG1115+080 suitable for 
time delay estimation; 

51) V band observations from November 1995 to June 1996 with four instruments: Hiltner 
2.4 m telescope, the Wisconsin-Indiana-Yale-NOAO 3.5 m telescope (WIYN), the Nordic 
Optical 2.5 m telescope (NOT), and the Du Pont 2.5 m telescope [48]. The light curves are 
presented in Fig. 1. 

52) R band observations with the SMARTS 1.3 m telescope at GTIO and the 2.4 m 
telescope at the MDM Observatory in 2004-2006 [29]. The light curves of S2 are presented 
in Fig. 2 with hollow circles. 

53) R band observations with the 1.5 m telescope at Maidanak observatory in 2001-2006 
[60]. The light curves of S3 are presented in Fig. 2 with filled circles. 

The time delay between Al and A2 spaced 0."48 apart is about several hours [6], therefore 
we average light curves of Al and A2 and compose the combined light curve A. 

The PG1115+080 variability in 1995-1996 (SI) is less than 0.2 mag, moreover authors 
themselves call their photometry data preliminary and note the presence of correlated errors, 
probably, resulting from the PSF determination error. In 2004-2006, the quasar variability 
amplitude was a bit larger then in 1995-1996 observation season and reached 0.4mag [61]. 
In the GTIO+MDM light curves (S2), the scatter in the data is too high, so the internal 
quasar variability becomes almost indistinguishable except Al and A2 components. Maida¬ 
nak observations (S3) are more complete and homogeneous data of long-term monitoring of 
the object. 
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Pic. 1. The light curves of PG1115+080 components in 1995-1996 based on Schechter’s data; the A 


light curve (averaged A1 and A2 light curves) is shifted by l.bmag 
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Pic. 2. The light curves of PG1115+080 components in 2001-2006 (hollow circles - CTIO+MDM 

data, filled circles - Maidanak data) 

To combine data from different instruments and in different bands, we do as follows; 

1) We convert all data from magnitudes to relative fluxes. It doesn’t affect the position 
of the maximum of the cross-correlation function. 

2) Then for each season we subtract average flux level from data to remove inter-season 
trends. It doesn’t affect the position of the peak of the cross-correlation function either, but 
essentially affects its width making it smaller and hence increasing accuracy of time delay 
determination, as linear trends and slow inter-season changes lead to high positive correlation 
in a wide range of delays (hundreds of days and more) and prevent us from obtaining time 
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delays with the required accuracy. Additionally, this alignment procedure allows exclusion 
of the trend caused by slow microlensing. 

3. STATISTICAL ANALYSIS OF THE PG 1115+080 LIGHT GURVES WITH THE 

MGGF TEGHNIQUE 

There are traditional cross-correlation techniques of the time delay estimation for evenly 
spaced time series. However astronomers usually have to deal with unevenly spaced time 
series. The problem of unevenness is solved in different methods in different ways. One can 
use linear or other interpolation and transform time series into evenly spaced ones and then 
use the classical procedure of the cross-correlation analysis. In doing so, the interpolation 
errors are introduced into initial data. Those errors can be mitigated if only one time set is 
interpolated leaving the other time set unchanged. A typical example of such an approach is 
the interpolation method of Gaskell and Sparke, hereafter called IGGF [14, 15]. In the discrete 
correlation function method of Edelson and Krolik , hereafter called DGF [28], the authors 
try to avoid interpolation, but they have to smooth and approximate the cross-correlation 
function itself. Another class of methods involves construction of some continuous curve, 
approximating quasar variability light curve, and htting it to both observational light curves 
through the minimization of the function to hnd parameters of that model curve and a 
time delay. 

The variety of techniques used for estimation of time delays in the PGl 115+080 and 
mentioned in the Introduction is due to the need to work at the limit of a method resolution. 
Our aim is to use available dehcient information as accurately as possible. 


3.1. The MCCF technique 


To analyze available light curves of PG1115+080, we use method that is intermediate 
between GGF and DGF and is based on linear interpolation of the light curves [15]. Let 
A{ti) and B{ti) be the measured fluxes of two components of a GLS at a time moment t*. 
The correlation coefficient between two time series A{t) and B{t) is calculated as follows: 


MCCF{t) 


J_ (Aj - A){Bj - B) 

M ^ 

*5 


O'aO'b 


( 1 ) 
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where M is the number of data pairs (Aj, Bj) for every r — Atmax < Atij < r + Atmax {Atij 
is a time lag between ti of the time series A and tj of the time series B), is a standard 
deviation, x is the mean of x. As the time series are unevenly spaced, we have to interpolate 
data at least in one of the light curves. If both time series have the same photometric accuracy 
and observations density, then initially we match the light curve A{t) with the interpolated 
light curve B{t + T) and calculate the cross-correlation coefficient for every r, next the light 
curves are interchanged, i.e. the observed light curve B{t) is matched with the interpolated 
light curve A{t + T). Thereby correlation coefficients are calculated twice for every time delay 
value, then average values of cross-correlation coefficients are calculated. 

We modihed the ICCF method by adding a restriction on usage of interpolated data: 
light curves are interpolated only at points which are sufficiently close to time of actual 
observations. The interval Atmax dehning such proximity is given a priori and is selected 
according to the characteristic time of the quasar internal variability and the need for trade¬ 
off between interpolation errors and a sufficient number of pairs of data points to calculate 
the cross-correlation coefficient. Unlike the ICCF technique, where the number of pairs of 
points in calculation of the cross-correlation coefficient is equal for all time shifts, in MCCF 
this number varies, as we neglect some of the data points. All points are used in the analysis 
only for the zero delay and when a time delay does not exceed Atmax- For time delays 
exceeding Atmax, fhe number of pairs are always less than the maximum one. Yet we need 
a sufficient number of pairs for the significant estimate of the cross-correlation coefficient. If 
one choose too low Atmax, the number of pairs of data points for calculation of the correlation 
coefficient for some time shifts may be insufficient. If the number of pairs for some time delays 
is less than ten (M < 10), we exclude the obtained value of the correlation coefficient from 
consideration. If one increase the interval Atmax, the number of pairs for analysis grows, but 
our method tends to the ICCF method, i.e. the contribution of interpolation errors grows. 

The ICCF technique implies not only interpolation, but also extrapolation; cross- 
correlated time series are extended at both edges with constant values equal to the first and 
last observation values. That means adding a large number of pairs of points not containing 
useful information, but smoothing the cross-correlation function, because inconclusive low- 
frequency part of the time series power spectrum is intensihed. We extrapolate light 
curves beyond observation seasons only within the interval Atmax- Those modihcations was 
suggested by Oknyanskiy [36], and the modified method was called MCCF (Modified method 



of Cross-Correlation Function). 

Light curves of GLQ components are irregular time series, but measurements of fluxes for 
all components are performed in the same days and in the same frames. Therefore correlated 
errors can enhance correlation near zero delay or weaken it if the GLQ components are closely 
spaced and there is a systematic anti-correlation. Some techniques, particularly the DGF one, 
propose to exclude data points obtained in the same night from the analysis. Unfortunately, 
we can not afford it, because then we have absolutely no data at the zero lag. Stability 
and robustness of the MGGF technique is repeatedly conhrmed: it has been successfully 
applied to determine the time delay between the components of the hrst gravitational lens 
QSO0957+561 [53], and to determine the time delay between flux changes in the near infrared 
and optical regions of the spectrum for several active galactic nuclei [35], and for measuring 
time lags between flux changes in the optical continuum and emission lines, near-IR and 
optical continuum for NGG 4151 [33, 34, 37]. It is also used to determine the delay in the 
gravitational lenses QSO2237+0305 (Einstein Gross) [24] and UM673 [25]. 

3.2. Application of the MCCF technique to light curves of PG1115+080 

Expected values of time delays in the gravitational lens PG1115+080 are of the order 
of dozens of days [30], that is signihcantly greater than the typical intervals between 
observations within a single season and less than an observation season duration, so chances 
of determining time delays, provided that there are an observable variability and lack of 
microlensing, are sufficiently large. 

To successfully calculate a time delay, intervals between observations should not exceed 
the characteristic time of the variability of the object and the expected time delay. For the 
considered light curves of PG1115+080 it holds true. An observation season should be long 
enough so that there were actual variations of the brightness of the quasar. The time delay 
determination can be complicated due to the variability caused by microlensing by the stars 
of the galaxy-lens. This variability has a characteristic time of a few months, and can be 
observed in the form of slow inter-seasonal trends and changes. In any case, inter-seasonal 
changes will only hinder the estimation of the time delay. 

In our case, when there are several observation seasons and gaps between seasons are 
signihcantly longer than the duration of the seasons, trends between seasons are not carry 
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any useful information and lead to high correlation coefficients over the entire range of time 
delays and low relative amplitude of the main maximum of interest. Therefore we align the 
average flux of the seasons to the same level to avoid inconclusive correlation. It can be done 
because our interval of interest is small enough and points from neighbour seasons do not 
overlap when we shift light curves relative to each other by a time delay. 

The cross-correlation functions, calculated with the MCCF technique for the light curves 
of B and C components from Maidanak observations in 2004-2006, are presented in Figs. 3 
and 4. The maxima of cross-correlation function for the seasons 2004, 2005 and 2004-2005 
coincide with each other within the conhdential interval, while the cross-correlation function 
for the 2006 has several peaks making it impossible to determine the maximum and therefore 
the time delay. 

Impact of data obtained in 2006 leads to the fact that the cross-correlation function for 
2004-2006 has no clearly dehned maximum. So if one determines the position of the maximum 
of the cross-correlation function using the median value, this position will be shifted towards 
smaller delays. It is likely that the impact of 2006 data resulted in tbc = 16.4 days, calculated 
by Vakulik team [61], that differs downwards from the values tbc = 23.7 and tbc = 25.0 
calculated by Schechter and Barkana [2, 48]. 

We noticed that the light curves of the Al, A2, B and C components in 2006 represent 
almost linear trend. And in the Al and C components, rapid oscillations are observed, which 
can be caused by observation errors. We believe that these data can reduce the statistical 
signihcance of estimated time delay values and even lead to biased results. Therefore we 
decided to exclude the 2006 season from further statistical analysis. 

Fig. 5 shows cross-correlation functions calculated for the B and C light curves in 1995- 
1996, 2004-2005, and the combined 1995-2005 light curves. It can be seen here that the 
position of the maximum of the cross-correlation function, and hence the delay tbc based on 
observations in 1995-1996, is greater than based on observations in 2004-2005. The combined 
1995-2005 light curves also yield the time delay value tbc = 22 days. This is the same value 
as based on observations in 2004-2005, as this observation season provides more pairs of 
points for the calculation of a cross-correlation coefficient. 

The cross-correlation functions calculated for pairs AB and AC in 2004-2005 are presented 
in Fig. 6. Positions of the peaks correspond to the delay of the internal quasar variability 
manifestation in A component relative to B component tba = 10 days and relative to C 
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Pic. 3. Cross-correlation functions for the BC pair of components based on S3 (Maidanak) 
observations in 2004, 2005, and 2004-2005. 
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Pic. 4. Cross-correlation functions for the BC pair of components based on the S3 (Maidanak) 

observations in 2006 and 2004-2006. 


component tac = 12 days. 

The CTIO + MDM data (S2) overlap with Maidanak data (S3), but the quasar internal 
variability is practically indistinguishable in light curves of B and C components due to the 
large scatter of the data points. Cross-correlation functions for the data shown in Fig. 7 do 
not allow determining the time delays. Vakulik and coauthors also reported that is was rather 
problematic to use the CTIO + MDM data to determine the time delays with their method. 
We have tried to include the light curve of the A component from S2 into the analysis in 
combination with Maidanak data, but it does not affect obtained time delay values. 
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Pic. 5. Cross-correlation functions for the BC pair of components based on the SI data (dotted 
line), S3 data (dashed line), combined SI and S3 data (continuous gray line) 
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Pic. 6 . Cross-correlation functions for the AB and AC pairs of the components based on S3 

(Maidanak) observations in 2004-2005. 

3.3. Uncertainties 

To calculate the uncertainty of time delay determination, we use Monte-Carlo approach 
that means application of the MCCF technique to a statistically large number of simulated 
light curve pairs with the specihed time delay. The simulated light curves have to reflect the 
statistical properties of the quasar variability and properties of the observed light curves. For 
example, one of the ways to assess robustness of the method is to simulate the light curves 
from the real data and compare the light curves with some points removed by chance [38]. 
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Pic. 7. Cross-correlation functions for the AB, AC, and BC pairs of the components based on S2 

(CTIO+MDM) observations in 2004-2006. 

However that approach underestimates time-delay uncertainties because not all statistical 
properties of quasar variability are incorporated into the simulated light curves. Therefore we 
decided to construct every pair of light curves based on a simulated light curve of the quasar 
variability with randomized parameters. The second light curve is produced by shifting the 
hrst one by the specihed time delay. Then we select points according to actual observations 
moments. 

The quasar PG1115+080 can be referred to the family of quasars with average statistical 
properties. A quasar power spectrum in optical band can be modelled as power law with a 
mean index a = —1 [67]. To simulate the quasar light curve, we need an algorithm that is 
able to produce time series with the power-density spectrum S{u) ~ u°‘ randomizing both 
the phase and the amplitude of the Fourier transform of a signal. Such an algorithm was 
proposed by Timmer& Konig [58]. 

The first step is to generate the power-density spectrum; 

^(cu) = iV(0 ,+ tNiO, l)./i (2) 

where A^(0,1) is a normally distributed stochastic variable with zero mean, and the standard 
deviation equal to 1. The simulated light curve is produced by applying the inverse Fourier 
transform to that power spectrum followed by renormalizing to the observed magnitudes. 
The second light curve is generated by shifting the simulated light curve by given time delay. 
To reflect the unevenness of actual light curves, the generated light curves are resampled 
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Pic. 8 . Histograms of time delays recovered from cross-correlation analysis of 1000 Monte Carlo 
simulations of the A and C, B and A, B and C light curves. 

as the observational ones. To simulate the photometry uncertainty, normally distributed 
random variables with zero mean and standard deviation equal to the actual photometry 
errors are added to the simulated light curves. 

The Monte-Carlo histograms are presented in fig. 8. We assume the range where 68 % of 
the 1000 realizations falls to be a confidence interval of the time delay determination (la). 

4. DISCUSSION 

Our results of time delay determination for PG1115+080 are summarized in the last row 
of the Table 1. The uncertainties correspond to la. Time delay values, calculated by other 
researchers, are also presented in the table for comparison. 

Our estimate of tba match the results of Schechter and Barkana within the accuracy of 
the analysis, but differs greatly from tba calculated by Vakulik et ah for the same Maidanak 
data except 2006 season that we excluded from the analysis. As for tac, almost all techniques 
yield similar time delay approximately 12 days. Only recent Rathna Kumar’ result is a bit 
larger. The time delay tbc, obtained by Vakulik et ah and Rathna Kumar et ah, is less than 
tbcj calculated by Schechter and Barkana and confirmed with our MCCF technique, by a 
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Table 1. Time delays between images of the lensed quasar PG1115+080 estimated with different 
techniques. 


References 

Season 

tba 

tac 

tbc 

tac 1 Tba 



(days) 

(days) 

(days) 


P. Schechter et al [48] 


14.3 ±3.4 

9.4 ±3.4 

23.7 ±3.4 

0.7 ±0.4 

R. Barkana [2] 

1995-1996 

11 . 7^|9 

13.310 

25.010 

1 1 q + 0.18 

-'-•■'-'^- 0.17 

J. Pelt et al [40] 


15.5 ±1.8 

10.3 ±1.9 

25.8 ±2.4 

0.7 ±0.2 

E. Eulaers & P. Magain [12] NMF 



11.7 ±2.2 

OQ 0 + 2.8 

^^•^- 3.0 


E. Eulaers & P. Magain [12] MD 



7.6 ±3.9 

17.9 ± 6.9 


V. Vakulik et al [61] 

2004-2006 

4.410 


16.410 

<^0 + 2.55 
^•*'^- 2.01 

S. Rathna Kumar et al. [47] 


8.9 ±8.4 

18.3 ±4.4 

13.2 ±9.0 

2.1 ± 2.5 

E. Shimanovskaya et al. [52] 

2004-2005 

loO 

I 2 I 4 

220 

1 9 + 0.4 
-'^•^- 0.5 


factor of 1.5. 

In the last column of the Table 1, the ratio tac/^ba for PG1115+080 is presented. The 
time delay ratio doesn’t depend on the Hubble constant and can be used for lens models 
comparison. Keeton and Kochanek investigated different models of the PG1115+080 lensing 
galaxy and theoretically calculated the ratio tacI^ba fo be in the range from 1.35 to 1.47 
[22]. As can be seen in the last column of the Table 1, the Schechter’s value tacI'Tba is 
signihcantly less than model predicted values, the ratio calculated by Barkana is close to 
the model predicted range, Vakulik’s value is more than twice as large as the theoretically 
predicted value and has the large uncertainty. Our value tac/tba = l-2^o;4e, is close to the 
model predicted value. 

Based on measured time delays in hve lenses (RXJ0911+0551, PG1115+080, 
SBS1520+530, B1600+434, and HE2149-2745), Kochanek estimated Hq = 

48^4/cm Mpc“^, if a lensing galaxy has isothermal mass distribution, and 
Hq = 71 ±6 kms~^ Mpc~^, if a lensing galaxy has constant mass-to-light ratio [23]. Recently, 
the values TTo = 69 ± ± km s~^ Mpc~^ [51] and Hq = 71.8 ± 5.8 km s~^ Mpc~^ 

[47] were derived based on lightcurves of 18 gravitationally lensed quasars. Those Hq 
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Table 2. Recent Hubble constant estimates. 


Method 

Ho, km s ^ Mpc ^ 

References 

Cepheids (HST) 

73.8 ± 2.4 (3.3%) 

Riess et al 2011 [45] 

CMB: WMAP+SDSS, 

70 ± 4 (5.7%) 

Tegmark et al 2004 [56] 

WMAP7, 

71.0 ± 2.5 (3.5%) 

Jarosik et al 2011 [19] 

PLANCK 

67.3 ± 1.2 

Planck Collaboration 2014 [41] 

Sunyaev-Zeldovich effect 

66 ± 11 (16.7%) 

Jones et al 2005 [20] 

AGN IR lags 

73 ±3 

Yoshii et al. 2014 [66] 

Gravitational lenses 

61 ±7 (11.5%) 

Courbin 2003 [9] 


68 ± 6 (8.8%) 

Oguri 2007 [31] 


69 ± 6^*“‘ ± (15%) 

Sereno & Paraficz 2014 [51] 


71.8 ± 5.8 (8%) 

Rathna Kumar et al. 2014 [47] 




values agree with a recent estimate of the Hubble constant based on Cepheids observations 
in the HST Key Project Hq = 73.8 ± 2Akm s~^ Mpc~^ [45]. Combining WMAP and 
SDSS CMB spectrum observations with the assumption of the flat Universe gives value 
Ho = 70t.o^tkms~^ Mpc~^ [56]. New results of the ESA’s Planck project testify to the 
smaller speed of expansion of the Universe Hq = 67.3 ± 1.2 km s~^ Mpc~^ [41], however in 
some recent works this result is called into question because of frequency channel calibration 
issues of the Planck telescope [54]. Another method of the Hubble constant estimation 

I 1 nrand \ Qsyst .. 

based on observations of Syunyaev-Zeldovich effect gives Hq = 59_grand L^syst km Mpc~ 
for standard cold dark matter with Qm = 1-0, Ha = 0.0 or Hq = 66'tio-s km s~^ Mpc~^ if 
Hm = 0.3, Ha = 0.7 [20]. Investigation of lens models shows that the Ho estimated based on 
gravitational lensing is close to values obtained with other methods only if a lensing galaxy 
contains small or zero portion of dark matter [18]. Hence, time delays between gravitational 
lens components allow estimating not only Hubble constant value, but portions of dark 
energy and dark matter also. 

The time delay tbc is less affected by systematic errors in comparison with shorter delays 
tba and tac- For this reason it is more often used for the Hubble constant estimation (e.g. see 
[22, 48]). Disregarding the whole variety of lens gravitation potential models, let’s take into 
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account just the fact that the time delay is inversely proportional to the Hubble constant 
[4, 8, 21, 43, 50]. The Hq value, estimated based on tbc calculated by Vakulik et ah or 
Rathna Kumar et ah, is 1.5-times greater [60] and hence closer to values measured with other 
methods, e.g. observations of Cepheids. However our value of tbc doesn’t support that result. 
Vise versa; it takes us back to the traditional situation when the Hubble constant estimates 
from time delays in gravitational lenses are signihcantly less than Hq values measured with 
other methods. Notice that the value of the Hubble constant measured in the Planck program 
is less than the Hq value derived from Cepheids observations as well [41]. 

5. CONCLUSION 

At hrst glance, the problem of measuring time delays between manifestations of internal 
quasar variability in components of a GLQ seems quite simple. However, actual observations 
are far from ideal case of evenly spaced time series. Except for observations uncertainties 
and large gaps in the light curves requiring usage of interpolation, superposition of micro- 
lensing and internal variability of the quasar and quite small amplitude of the quasar internal 
variability pose many problems. Even for such ’convenient’ lens for time delay estimation as 
PG1115+080, different researchers obtain different values of time delays, as indicated in the 
Table 1. This is particularly due to work at the limit of resolution of time delay methods, lack 
of observational data needed for robust estimation of time delays, and some parts of data 
not containing useful information for time delay estimation, as in the case of 2006 season. 

Different values of time delays lead to dissimilarity in the Hubble parameter estimation. 
Values of this universal constant calculated from time delays between GLS components are 
usually smaller than values obtained with other techniques. New estimates of time delays in 
the PG1115+080 gravitational lens, published by Vakulik et al, suggest larger value of the 
Hubble constant then previous estimates [60] narrowing a gap between values obtained with 
different approaches. 

We analysed all available optical observations of the gravitational lens PG1115+080 
with our cross-correlation technique to conhrm or argue against that result. We failed to 
determine time delays based on combined data because of inhomogeneity of observations 
and possible errors of their initial processing. Observations of PG1115+080 in R band on 
1.3-m SMARTS telescope (GTIO) and on 2.4-m telescope of MDM Observatory in 2004- 
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2006 [29] are unsuitable for the time delay estimation with both Vakulik’s [61] and ours 
techniques, because the quasar internal variability is bleared by chaotic abrupt behavior 
of light curves that may be caused by microlensing and photometry errors. Among three 
available long-term data sets of optical observations of PG1115-I-080, the most complete 
and homogeneous ones were obtained at Maidanak observatory (Uzbekistan) in 2004-2006 
[60]. But we hnd that light curves of all four components of the GLQ in 2006 represents 
almost linear trends. Fast variations are observed only in Al and G components that can be 
as due to photometry errors or microlensing. These fast variations can decrease statistical 
signihcance of the time delay estimates or even produce misleading results. Therefore we 
decided to exclude observations in 2006 from the cross-correlation analysis. Time delays 
estimated with the MGGF technique [36] applied to Maidanak observations in 2004-2005 
are in better agreement with previous results by Schechter [48] and Barkana [2] obtained 
for 1995-1996 light curves with two different techniques of statistical analysis than with new 
estimates by Vakulik et al [60, 61] and Rathna Kumar et al [46]. The ratio tag/^ba is in 
agreement with Barkana’s value and model predictions supporting adequacy of our results. 
Our estimates of time delay suggest smaller value of the Hubble parameter than Vakulik et 
al results. 

As long as different methods yield different values of time delays for the same data, it’s 
early to talk about estimation of the Hubble constant. We have to identify reasons of such 
difference that can be due to interpolation and accounting for microlensing. To attain better 
accuracy, new homogeneous long-term observations and new methods of analysis taking into 
account microlensing and weight of every data point are also needed. 
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